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O ■ ABSTRACT 

Remnants of AGN jets and their surrounding cocoons leave colossal magnetohydro- 
dynamic (MHD) fossil structures storing total energies ~10 60 erg. The original active 
galacic nucleus (AGN) may be dead but the fossil will retain its stable configuration 
resembling the reversed-field pinch (RFP) encountered in laboratory MHD cxperi- 
\ ments. 
Q_c In an RFP the longitudinal magnetic field changes direction at a critical distance 

from the axis, leading to magnetic re-connection there, and to slow decay of the large- 
scale RFP field. We show that this field decay induces large-scale electric fields which 
can accelerate cosmic rays with an E~ 2 power-law up to ultra-high energies with a cut- 
\ off depending on the fossil parameters. The cut-off is expected to be rigidity dependent, 

implying the observed composition would change from light to heavy close to the 
cut-off if one or two nearby AGN fossils dominate. Given that several percent of the 
universe's volume may house such slowly decaying structures, these fossils may even re- 
energize ultra- high energy cosmic rays from distant/old sources, offsetting the "GZK- 
losses" due to interactions with photons of the cosmic microwave background radiation 
and giving evidence of otherwise undetectable fossils. In this case the composition 
■^j- ■ would remain light to the highest energies if distant sources or fossils dominated, but 

otherwise would be mixed. It is hoped the new generation of cosmic ray experiments 
such as the Pierre Auger Observatory and ultra-high energy neutrino telescopes such 
as ANITA and lunar Cherenkov experiments will clarify this. 
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1 INTRODUCTION the highest energies are, however, subject to interaction 

with the cosmic microwave background radiation (CMBR). 
The energy spectrum of cosmic rays (OR) extends from be- . , , r , , r , , 

, . , _,„ 2 n ,r m, i n r i ■ 111 the rest frame of a UHE OR proton, photons of the 

low 1 Gev up to at least 10 ev. the bulk oi the cosmic „ , r „, mn , , , , ... , . 

i. 2.73 K CMBR are strongly blue-shifted to gamma-ray 

rays below ~10 ev are probably Galactic. It is currently . „, ., , „ . „,,„.,, . , ,. 

energies, the thresholds tor Bethe-Heitler pair production 
unknown whether the ultra-high energy (UHE) cosmic rays , . , . , ,. , ITITn " 

° , , , an( i pion photo- production by UHE OR protons on the 

are Galactic or extragalactic in origin and whether they are CMBR are close to 2 x 10 17 eV and 2 x 10 19 eV, respectively, 
accelerated or result from emission by topological defects protons at 3xl() i9 eV and 3xl0 20 eV t icall lose a } 

or decay of super-massive particles. If they are acceferated, . ,. r . . .. r -, ^ , , n , „n s , 

, , fraction of their energy in a time of 1 Gpc/c (3 x 10 y) and 

they are almost certainly extragalactic in origin. This is be- 1r , ,, , ,„ ,„7 ^ , . , „, r . 

10 Mpc/c (3 x 10 y), respectively, the importance of pion 

cause their gyroradn are sufficiently farge that an amsotropy , , , , , „, T „„ n . , , , ^ . 

. , photoproduction on the CMBR was hrst noted by Greisen 

would be expected toward a galactic source population (not (1Q66) and Zatsepin & Kuzmin (lg66) and ^ 

observed), and because there are no obvious Gafactic source ,. , , . r , . mi mi 

. . ' „ r they predicted is referred to as the GZK cut-off . the 

popufations having magnetic helds such that gyroradn of . , , . ,. . , ,. . n „ .,.„ 

f TT ^_, pion photoproduction interactions resuft m a flux of UHB 

UHE CR are compatible with their size. , . ,. , , , it^ms j. • » i • i 

neutrinos, often referred to as GZK neutrinos , which 
Possibfe extragalactic acceleration sites include , , ,i r , rrn 

would point back to halos surrounding the UHE cosmic 
hotspots of giant radio galaxies, the intergalactic medium, „ , . , , , ,„„„„, 

° ray sources. Several experiments (lakeda et al. (2003), 

gamma ray bursts and blazar jets. Extragalactic CR at , . n r ^ nn r\ ^ n i Air? i 

& ■' j a g lrd et a i_ (1995J ; Connolly et al. (2006), Abu-Zayyad et 

al. (2005)) have reported UHE CR events with energies 

. _ ., „ _ r .» . , above 10 20 eV. The Auger array (Auger Collaboration 
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t E-mail: rprother@physics.adelaide.edu.au (2004)) will soon report its results with excellent statistics. 
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Meanwhile, experiments such as ANITA (Barwick et al. 
(2006)) and lunar Cherenkov experiments are in progress 
to detect the UHE neutrinos which should give clues to 
the origin of UHECR. See Protheroe and Clay (2004) and 
Falcke, Gorham & Protheroe (2004) for recent reviews of 
UHECR and radio techniques in UHE neutrino and particle 
astronomy. 

Remnants of jets and their surrounding cocoons may 
persist long after their parent AGN fade from view. These 
colossal MHD structures decay slowly and yet may retain 
their relatively stable self-organized configurations. Decay 
depends on the structure circuit resistance, and lifetimes 
could be quite long, given the large inductance of the cir- 
cuit, an initial outward current along the jet and a return 
current back along an outer sheath or cocoon around the jet. 
The current closure condition requires that a large fraction 
of magnetic fields resistively dissipate (Lesch et al. 1989). 
Also Lyutikov & Ouyed (2007) argue for direct inductive 
acceleration in sheared relativistic jets while the jet are on, 
though not afterwards. We concentrate here on the fossil era, 
when a substantial fraction of jet energy has already been 
invested in the magnetic energy of the fossil and induction 
occurs through decay, not in the earlier building period - 
though both eras can contribute to acceleration schemes. 

Helical structures are common in AGN jets, arising dur- 
ing jet formation from rotation of magnetized plasma accret- 
ing toward the central black hole - this could be enhanced 
in the case of binary black hole systems. Azimuthal electric 
currents are therefore likely, yielding a magnetic field com- 
ponent along the jet direction. In radio galaxy jets boundary 
layers are often clearly visible in radio polarization (see, e.g. 
Laing & Bridle (2002) for observations and interpretation of 
the jets of 3C 31). Ostrowski (2002) has argued that bound- 
ary layers in jets are sites of particle acceleration. 

Laboratory MHD experiments show that reversed-field 
pinches are fairly stable. On the immense scale of these "fos- 
sil jets", the decay time from instability can be billions of 
years. However, decay of these colossal MHD structures on 
such time-scales result in electric fields capable of acceler- 
ating existing populations of lower energy cosmic rays up 
to ultra high energies with a flat spectrum extending from 
some minimum rigidity (momentum/charge) determined by 
fossil dimensions, magnetic field and decay time. We first re- 
view the properties expected of large magnetic fossils, based 
on extensive work in the plasma physics community over the 
last four decades. 



2 STABILITY OF MHD STRUCTURES 

We assume that the governing, evolved equilibrium of ini- 
tially jet-driven, then later, long-lived magnetic structures 
will be a reversed field pinch (Taylor (1986), Bellan (2000)). 
This axisymmetric configuration has magnetic field compo- 
nents (in cylindrical coordinates - r, 4>, z) 

B z ~ B ~ y^P (1) 

in an MHD equilibrium made stable by optimally efficient 
radial profiles, demanding a minimum of B z - B^, is built 
by the jet current. 

In a reversed-field pinch the crucial geometric feature is 
the longitudinal magnetic field, which changes direction at 



some critical distance from the jet axis. In experiments the 
RFP can be curved into a torus, creating a spheromak (Bel- 
lan (2000)). We use a simple cylindrical RFP model below 
because, though a nearly spherical geometry may be better 
for an evolved fossil (especially in rich clusters with nearly 
isotropic pressures) , the spheromak geometry imposes many 
field curvatures effects in particle orbits which are a compli- 
cation beyond the essential physical point of the reversed 
axial field. 

Stable, low /?(/? = 2fi P/B 2 where P is the plasma 
pressure) reversed field pinch radial profiles have several 
properties: 

(i) B z reverses near the outside of the confinement region. 
This is the crucial shear that stabilizes modes that can inter- 
change flux tubes of plasma. It also prevents kink (m = 1), 
current-driven instabilities. 

(ii) To suppress "sausage" modes driven by pressure, a 
value of /3 < | is essential. The smaller the /3, the more 
stabilized a confinement geometry becomes. 

(iii) A conducting "wall" close enough to the plasma core 
to suppress the kink, m = 1, current-driven internal kink 
modes. This also completely damps all external kink modes. 

(iv) A pressure profile P(r) that is hollow or very flat in 
r, to suppress interchange of magnetic flux tubes near the 
magnetic axis. 

Melrose, Nicholls & Broderick (1994) showed that com- 
plicated equilibrium models of an isolated, force-free mag- 
netic flux tube could be constructed with cylindrically sym- 
metric, force-free magnetic fields. These had a sequence of 
concentric layers with piecewise-constant a, where a is a 
running index of the equilibria in concentric cylinders, each 
obeying the force-free condition, j(r,<f>,z) — a(r)B(r,<f>, z). 
They used analytical methods to construct models of the 
surface currents on these equilibria, finding some conditions 
such that the net current flux in the tube was zero. 

For example, a "two-a model" permits current neutral- 
ization with a non-zero magnetic flux. An infinite series of 
equilibria are possible, yielding stable reversed-field geome- 
tries with various boundary conditions. Often a nearby con- 
ducting boundary simplifies the equilibrium as an ideal ap- 
proximation. 

Are such equilibria at energy equipartition between 
magnetic and relativistic plasma energies? Dunn & Fabian 
(2004) studied in cluster "ghosts" the limits on k/ f, where k 
is the ratio of the total relativistic particle energy to that in 
electrons radiating between 10 MHz to 10 GHz, and / is the 
volume-filling factor of the relativistic plasma. Strikingly, 
none of their "bubbles" had a simple equipartition between 
the pressures from the relativistic particles and the mag- 
netic field. Furthermore, k/f did not strongly depend on any 
physical parameter of the host cluster. Though at first there 
seemed to be two populations - k/f values around 2, another 
bunch around 300 - this did not hold up (Dunn, Fabian & 
Taylor (2005)). The apparent bimodality of the k/f distri- 
bution could be explained as arising from two kinds of jets 
- electron-positron, giving a low value for k, and electron- 
proton. If protons are the extra particles needed to maintain 
pressure equilibrium, but unseen in the radio emission, k is 
high. Also, bimodality could come from either a non-uniform 
magnetic field, or a filamentary structure in the lobes. 

Both these possibilities are consistent with a reversed- 
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field equilibrium, since fields vary spatially. Variations in re- 
acceleration, which may occur from reconnection events in 
the field-reversal zone, could also affect the k/f measure. At 
this point we know too little to infer much. Clearly, though, 
constraints on the magnetic field of cluster bubbles found by 
comparing the synchrotron cooling time to the bubble age 
show that no bubbles in the sample are in equipartition. 

This is plausible as laboratory experience with long- 
lived "magnetic bottles" such as tokomaks and reversed field 
pinches show that only those far from equipartition survive 
long. Only when the plasma pressure is far less by several 
orders of magnitude than magnetic field energy density do 
structures long persist. This is an important lesson of the 
applied fusion research program, now over half a century 
old: only low (3 is usefully stable. 

In this spirit, we shall use the simplest possible current- 
neutralized, reversed-field model, with a single alpha and 
a nearby ideal conducting wall. It has zero total azimuthal 
magnetic flux since some current is carried in our assumed 
conducting wall. This simplification allows us to concentrate 
on the particle acceleration mechanisms in a simple geome- 
try. Without this simplification, we would need a more com- 
plex equilibrium. 

2.1 Evolution of Reversed-Field Magnetic 
Structures 

Recent detections of several ghost cavities in galaxy clusters 
(Clarke et al. (2006)) - often, but not always, radio-emitting 
- suggest that the cluster hot plasma stays well separated 
from the bulk of the relativistic plasma on a timescale of 
~100 Myr. Some leakage is plausible, but robust survival 
does occur. In cases where the radio-emitting cosmic-ray 
clouds are confined within the observable hot gas reservoir of 
the cluster, the dimensions of the X-ray hole combined with 
the calculated gas pressure gives a "laboratory-like" mea- 
surement of the central black hole's energy that is converted 
into PdV work against the intra-cluster medium (ICM). An 
example is the X-ray cluster MS 0735.6+7421 observed in 
X-rays (Chandra) by McNamara et al. (2005) in which two 
large ~200 kpc cavities are seem in the X-ray emission; ra- 
dio emission comes only from the cavities. The total energy 
from the AGN outburst that created these giant magnetic 
bubbles, and probably re-heated the cluster gas, is estimated 
to be ~ 6 x 10 61 erg, of which a substantial fraction must 
be in magnetic field based on pressure arguments. Dynam- 
ical changes will drive inductive electric fields (Ensslin et 
al. (1997), Kronberg et al. (2004)). Referring to galaxy clus- 
ter embedded sources, Kronberg (2005) states that pressure 
balance between the magnetic field and gas (thermal and 
CR gas) within the holes, and the external thermal hot gas 
pressure implies typical B values of « 20-30 fiG within the 
holes. 

So, these ghosts may be RFPs. How have they survived 
so long? Benford (2006) proved a theorem: that for such jet- 
built RFP structures, the same simple MHD stability con- 
ditions guarantee stability, even after the jet turns off. This 
means that magnetic structures made stable while a jet is 
on can evolve into fossils that persist long after the build- 
ing jet current has died away. These may be the relic ra- 
dio "fossils" , "ghost bubbles" or "magnetic balloons" found 
in clusters and made visible by contrast against the X-ray 



emission as seen in Hydra A by Wise et al. (2007). Giant 
radio galaxies such as Cyg A are building such structures 
now. Such fossils have a massive inventory of magnetic en- 
ergy that can be ~10 60 erg (Kronberg et al. (2004)). Such 
enormous energy can only come from the gravitational infall 
energy of a supermassive black hole, when jets convey a few 
percent of the energy outward. 

These equilibria will evolve under the magnetic tension 
that sustains them. They are active, responsive structures, 
not a disordered ball of magnetic fields. In rich clusters these 
balloons may contract under external gas pressure. They can 
later rise in the gravitational potential and then expand in 
the lesser external pressure farther from the cluster core. 
Volume change will alter the plasma pressure as p dV work 
changes the magnetic equilibrium. This in turn drives in- 
ductive electric fields, which can accelerate ions. We envi- 
sion such inductive changes in response to external forces or 
internal evolution as the deep cause of high energy particle 
acceleration. 

The best examples of radio bubbles are 3C 84, M87 
and Her A. The edges of these bubbles are amazingly sharp, 
even in projection, as would be expected at the edge of a 
long-lived self ordered magnetic structure. Indeed, because 
of the helicity present in the magnetic field of the original 
AGN jet, the magnetic bubble it produces will retain this 
helicity and evolve to a self-organized force-free structure 
such as an RFP or sheromak, which are hard to pull apart 
- the toroidal field lines opposing the the lateral exansion of 
poloidal. The magnetic fields of these structures will there- 
fore undergo minimal mixing with the overall cluster mag- 
netic field, giving rise to the sharp edges seen radio images. 
They should thus also have larger magnetic fields than the 
surrounding intra-cluster medium. 

Initially stable reversed field pinches become endan- 
gered by the current- driven global kink instability. This 
seems to be the most likely way for structures to fail when 
the jet current source is on, since there is more free energy 
to drive instability. If a structure survives the current-driven 
era, later stability seems more probable, as loss of the jet 
removes free energy. But suppose the system fails to shed 
magnetic helicity K as dissipation of magnetic energy pro- 
ceeds (by reconnection, primarily). Such helicity loss is well 
known to lead to kink instability. Here we can turn to abun- 
dant laboratory experience as a guide. 

2.2 The Laboratory to Astrophysical Analogy 

Laboratory experiments provide much of our lore about re- 
versed field pinches. These are all low pressure devices (low 
/?). Taylor (1963) argued that in all systems with (3 < 1 mag- 
netic reconnection conserves global magnetic helicity, allow- 
ing stable geometries to evolve without losing stability (see 
also Taylor (1986)). This is so even without the infinite con- 
ductivity assumption that is common in astrophysics. "Un- 
freezing" the magnetic field lines from plasma demands re- 
sistivity and thus some form of magnetic energy dissipation. 
Even so, helicity can be preserved, as in the case of dissi- 
pation by reconnection. Thus /3-values approaching 1 may 
occur in long-lived systems if helicity (twist) may be shed 
through turbulent processes other than reconnection. Even 
if the plasma (3 in the jet is not small, after the jet switches 
off and develops into as fossil jet with a self-organized mag- 
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netic structure, plasma cooling and expansion is expected 
to lead to lower /3-values. Recently, Lapenta & Kronberg 
(2005) have compared detailed MHD simulations with radio 
observations of the jet of 3C 303, a FR II with a knotty 
structure, with the knots appearing like self-contained mag- 
netic structures launched along the jet. They find that in 
the knots (i ~ 10~ 5 , and that the jet is almost certainly ele- 
cromagnctically dominated rather than plasma dominated. 

A magnetically driven helical dynamo can also be 
described as "dynamical magnetic relaxation". Magnetic 
relaxation describes how structures in magnetic-pressure- 
dominated environments evolve to their equilibrium states. 
The fully relaxed state is the Taylor state, determined by 
minimizing the magnetic energy subject to the constraint 
that magnetic helicity is conserved. This leaves a force-free 
helical configuration with the field relaxing to the largest 
scale available, subject to boundary conditions. However, 
Taylor's theory by itself is not a dynamical theory since it 
does not provide a time-dependent description of how the 
large scale magnetic helicity evolves. 

We can draw an analogy between spheromak or 
reversed-field pinch formation and coronal field relaxation: 
the laboratory configurations form because the injection of 
magnetic helicity into a system takes the system away from 
the relaxed state. But helicity injection also drives tearing 
mode fluctuations that yield a finite, turbulent electromo- 
tive force. This in turn drives the system to a new relaxed 
state. In astrophysical coronae the situation is similar. The 
helical field produced by a velocity-driven dynamo inside a 
rotator gets buoyantly fed from below into the corona. This 
takes the magnetically dominated corona away from the re- 
laxed state, but fluctuations arise, driving the system back 
to a relaxed state. The corona thus acts like a "laboratory" 
plasma volume injected with helicity from below. Injection 
is not at one site, but throughout the coronal-surface bound- 
ary. While each of the individual injection sites (flux tubes) 
is like a single laboratory configuration if the field twists 
get injected along the tube, each tube can also open up to 
even larger scales by interacting with other tubes. Magnetic 
energy then appears on much larger scales - a difference 
between the laboratory and astrophysical cases. 

The dynamo effect in laboratory plasmas is a self- 
organization process. Resistive diffusion evolves the plasma 
away from a preferred state. The dynamo forces the plasma 
back toward the preferred, self-organized state. The result- 
ing magnetic configuration has excess free energy leading to 
instabilities (or turbulence) that relax the plasma toward 
a state of lower magnetic energy. The plasma relaxes to 
the lower energy state through the "dynamo" effect. The 
self-generation of plasma current reconfigures the large-scale 
magnetic field to one with lower energy. 

The structure of the relaxed state is partly captured 
by minimizing the magnetic energy in the plasma volume 
(Ub = J (B 2 /2fio)dV) subject to the constraint of constant 
magnetic helicity (K = J A-B dV), where A is the magnetic 
vector potential). Magnetic helicity topologically measures 
the knottedness of the magnetic field lines. The minimiza- 
tion yields a magnetic field such that the ratio of the current 
density to magnetic field, j/B, is a local spatial constant (the 
Taylor state; Taylor (1963)). For a recent comprehensive dis- 
cussion of these effects, which subsumes much of the earlier 



discussion especially as it unites theory and experiment, see 
Bellan (2000). 

An experimental test of the Taylor conjecture inferred 
the change in magnetic energy and helicity during a relax- 
ation event (Ji, Prager & Sarff (1995)). An approximate 
measurement modeled the instantaneous plasma state as a 
slowly varying MHD equilibrium. Since the fields during re- 
laxation change slower than an Alfven time, the fields satisfy 
the MHD force balance equation, 

VP = jxB (2) 

through discrete field twisting and braiding action. Solution 
of this equation with experimental constraints found that 
during relaxation the magnetic energy reduces by about 8%, 
while the magnetic helicity drops by about 3%. 

The Taylor conjecture gives a useful, experimentally 
verified framework to depict approximately the final state 
of the relaxation, or dynamo, process. But it provides no 
information on the physical mechanism of the dynamo. So- 
lution of the resistive nonlinear MHD equations reveals the 
detailed dynamics. The spatial fluctuations that underlie 
the dynamo are tearing instabilities driven by spatial gradi- 
ents in j/B, making the field lines tear and reconnect. The 
MHD description of such spontaneous magnetic reconnec- 
tion has been investigated for several decades. Linear theory 
shows that tearing instabilities grow on a timescale between 
a short Alfven time (L/va), and the long resistive diffusion 
time (noL 2 /ri, where r\ is the resistivity). For experimental 
plasma parameters, the theoretical growth time is approxi- 
mately 100 /xs, comparable to the observed relaxation time. 
In astrophysical jet-built structures, this can be billions of 
years. 

Tearing arises from a resonance between the fluctuat- 
ing and mean magnetic fields. In laboratory plasma tearing 
occurs at many radii within the plasma, leading to large- 
scale reorganization. For stationary and homogeneous tur- 
bulence, magnetic helicity K is relatively conserved com- 
pared to magnetic energy during relaxation. These dynamo 
effects play important roles, conserving the total magnetic 
helicity, except for resistive effects and a small battery effect. 

In astrophysical magnetic structure relaxation, charge 
and current imbalances adjust through the bulk, low energy 
plasma, which is collisional. The higher energy population 
can take advantage of the electric fields generated by decay, 
for they flow unimpeded along field lines, driven steadily by 
the electric fields. This is the essential picture in our model. 

The turbulent dynamo effect converts magnetic helicity 
from the turbulent magnetic field to the mean field, when the 
turbulence is electromagnetic. When the turbulence is elec- 
trostatic or due to the electron diamagnetic effect, the mag- 
netic helicity of the mean magnetic field gets transported 
across space. Two types of electric fields get generated dur- 
ing reconnection in the RFP. One is the helical electric field 
accompanying the helical tearing mode that causes the re- 
connection. The second is a mean, or axisymmetric, electric 
field induced by the flux change driven by the tearing mode 
dynamo. Clearly this is a complicated phenomenon. 

We do not attempt to treat such details. Instead, we 
shall simple take the inductive electric fields to be a frac- 
tion of the mean magnetic field. In the absence of detailed 
knowledge of the structure, this seems a good way to begin. 
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Figure 1. Magnetic field components of an RFP. 

2.3 Properties of a simple reversed field pinch 

The simplest idealization would be for an infinite cylindri- 
cal jet where the magnetic field in cylindrical coordinates 
(r,<t>,z) is (Lundquist (1951)) 

B r = 0, B4r) = B Ji(ar), B z (r) = B J (ar), (3) 

and this has been shown (Voslamber & Callebaut (1962)) to 
be stable for ar < 3.176. This radius for stability which we 
take to be R — 3.176/q is where a conducting wall with a 
large inertial mass would be present in an experimental sit- 
uation, and provide part of the circuit along which a return 
current could flow, which we assume here to be the cocoon. 
The magnetic field is shown in Fig. 1. Notice the longitu- 
dinal field changes sign at r cr i t = 2.405/a, the first zero of 
Jo(ar). 

One can use Ampere's law to find the current density 



-, s 1 B(r)a , , _o. 

j(r) = ( V xS) = -^- (Am 2 ), 



(4) 



which is everywhere proportional to the magnetic field. From 
this, the vector potential is 



A(r,4>,z) = -B(r,<p,z). 



(5) 



The outward electric current (positive j z ) is associated with 
radii r < r cr it and the return current with radii r > r cr i t , 
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(6) 
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and so some of the return current must be at radii r > R — 
in the outer part of the cocoon or in the interstellar medium 
of the host galaxy in order to exactly balance the outward 
current. 

The plasma pressure may be calculated from Eq. 1 of 
Robinson (1971), i.e. 



dP dB z B* d 
ar ar r ar 



(8) 



Solving this for B given by Eq. 3 results in P=constant. 
That is, for B given by Eq. 3 the tension of the azimuthal 
field lines perfectly balances the radially outward pressure 
of the longitudinal field. 



3 CHARGED PARTICLE MOTION AND 
ACCELERATION 

The critical radius, r cr it, where the longitudinal field B z 
changes sign, is also the radius where the longitudinal cur- 
rent density j z changes sign. It is also here that reconnec- 
tion may occur and may be a site of particle acceleration. 
Lesch, Birk & Schopper (2002) have considered acceleration 
in AGN jets due to a reconnection electric field along the 
length of the jet of an active quasar. That paper shows ac- 
celeration of hadrons to UHE should be possible in active 
AGN jets. The length of the jet is critical in their model, 
determining the maximum energy. We differ in that the re- 
connection electric field is azimuthal, i.e. around the fossil 
cylinder. Also, there should be many nearby remnants, so 
they can accelerate UHECRs, whereas active quasars tend 
to be at high redshift, and so their UHECRs will suffer the 
GZK cutoff if they are more than ~10 Mpc away. 



3.1 Energy of Reconnection: Plasma Wave Model 

Rather than assign a fixed fraction of the ambient magnetic 
field energy to the generated reconnection electric field, we 
can use results from laboratory observations in reversed field 
pinches (Ji et al. (2004); Ji, Prager & Sarff (1995)). The 
zone where B z reverses will probably be the reconnection 
site, since there only the azimuthal field needs to reconnect. 

In reconnection regions, micro-instabilities enhance re- 
sistivity, driving the reconnection rate. These arise from 
both electromagnetic and electrostatic modes. The electro- 
static lower hybrid drift instability has been invoked, but re- 
cent work shows that it does not occur in the high-/3 regions 
where reconnection occurs. There, electromagnetic instabili- 
ties prevail, as the B z field moves radially toward the recon- 
nection region, in our assumed cylindrical geometry, pref- 
erentially annihilating in the zone where B z reverses sign. 
This drives an azimuthal electric field. The plasma is highly 
collisionless, so the induced resistivity greatly exceeds the 
Spitzer level. The cause of this is ultimately the enhanced 
modes, such as the modified two-stream instability, which in 
experiment appear localized in the high-beta region where 
reconnection occurs. Magnetic fluctuations B\ appear also, 
correlated strongly with resistivity enhancement (Ji et al. 
(2004)). 

We can estimate the equivalent electric field generated 
in reconnection by noting that a momentum p per unit vol- 
ume appears in the electromagnetic waves, 



p/V = 



k{Bj) 



(9) 



Momentum moves from background plasma electrons 
and ions to the electromagnetic fluctuations at a rate 27, 
where 7 is the nonlinear damping rate on the electron-ion 
background plasma. In this highly nonlinear state, experi- 
ments show that 7 ~ u), the unstable mode frequency. Loss 
of momentum is a force equivalent to an effective electric 
field in the reconnection region, E ICC , given by 

(Bl) 



enE rcc = 2kj- 



W/Llo 



(10) 



The magnetic wave fluctuation levels seen in experi- 
ment are of magnitude B\ ~ 0.05 the ambient averaged 
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inflow field levels, B z . We thus set B\ = gB z and the 
fraction g — 0.05. Further, from experiment the observed 
wavelengths of this turbulence are the size of the coherence 
lengths in the fields, strongly suggesting a cooperative role 
(Ji et al. (2004)). The relevant wavenumber k is the inverse 
Larmor radius of background electrons in the plasma, 



tl = (10 m) x (T cV r 1/2 { ^ 



Bo 



(11) 



(Ji et al. (2004)). 

Scaling the results of laboratory experiments, we find 

E ICC = (1.5 x 10" 6 V m- 1 ) x * x 

'(■mf^iw^)' 1 < 121 

with n being the plasma density, and taking B z ~ 10 /iG and 
assuming the heated region has keV energies. To attain fields 
of order 10~ 4 V/m then requires g a bit larger, and perhaps 
B of 100 /xG. Plasma temperature may also be higher than 
keV, and density lower, also increasing E ICC . These parame- 
ter ranges would yield 10 20 eV protons after several encoun- 
ters with such structures of scale ~Mpc. Thus, we consider 
the possibility of UHE CR getting sequential boosts in en- 
ergy as a result of multiple encounters with different fossils 
in a similar way that Ip & Axford (1992) envisage Galactic 
cosmic rays get re-energized by multiple interactions with 
supernova remnants in the Galactic disk. 

We expect it to be realistic to assume that electric fields 
outside the reconnection zone can also accelerate particles. 
The cold plasma (pressure is low) responsible for currents 
which maintain the magnetic structure cannot short out 
these electric fields, since they are inductively driven ev- 
erywhere in the structure, allowing acceleration outside the 
reconnection zone. Next we investigate whether or not this 
alone will be viable for accelerating UHECR. In what fol- 
lows we assume that electric fields outside the reconnection 
zone can accelerate particles. 



3.2 Global Electric Fields from RFP Decay 

Electric fields from reconnection are emjs induced accord- 
ing to Faraday's law, and so the electric field will be more 
extensive. Assuming a flow of flux lines toward the recon- 
nection region, the field will be changing everywhere and 
will induce an electric field. The simplest way of estimat- 
ing this global reconnection electric field is by assuming an 
exponential decay of the magnetic field, 



B{t) = B o e~* /td< = c [Ji(ax)0 + J (ar)z]. 
Then, 

B 



_dA 
dt 

3.18 x 10 



B 



(13) 



(14) 



R 



( tdec \ 



10 ,uG J \ W0 kpc J \Gyr / 
x [Ji(ar)4>+ Jo(ar)i]e"* / * doc (V/m) (15) 
The rate of change of the particle momentum p — (3^mc 
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Figure 2. Three orthogonal views showing a typical orbit in the 
RFP fields, and critical radius (dashed). 



(neglecting losses) of a relativistic charge Ze is then 

$ = ZecB0x B + f3 dcc B), where /3 dcc = —^,(16) 

at Ctdec 

such that a simple 1st order integration scheme might be 



r(t + At) = r(t) +cAt/3(t), 
cAt 



(17) 



p(t + At) 



where r gyro (r) 



m + 



■p x 



r gyro(^) 

X0(t) x B(r) + f3 dcc B(r)], (18) 
Ifimc (lg) 



ZeB(f) ' 

Using this scheme we have conducted simulations of charged 
particle motion in the RFP magnetic field including the ef- 
fect of energy change in the induced electric field. We in- 
ject particles uniformly and isotropically over the surfaces 
of disks of radius R at both ends of the jet, and follow their 
motion until they escape. A typical example is shown in 
Fig. 2. 

Ultra-relativistic particles of charge Ze are injected with 
energy Eo (note for ultra-relativistic particles E « pc) and 
their final energies are binned as shown in Fig. 3. Since the 
induced electric field is in the same direction as the magnetic 
field, energy is gained as particle move along field lines. Thus 
we have also carried out simulations in which the particle's 
position is simply advanced along a field line, gaining mo- 
mentum ZeEAt in step cAt. Results for this case with high 
statistics (solid curves) are in agreement with the results for 
full gyro motion. We have shown that for a decaying RFP, in 
which the electric field is everywhere parallel to the magnetic 
field, cosmic ray nuclei will gain energy unimpeded except 
by pion photoproduction on the CMBR - as they undergo 
helical motion around magnetic field lines their momentum 
component parallel to B increases, while their momentum 
component perpendicular to B and the gyroradius remains 
constant, as is seen in the trajectory shown in Fig. 2. Note 
that the "double peaked" structure (Fig. 3) is due to sepa- 
rate contributions from injection at r < r cr it and r > r cr it, 
as we shall show in the next section. 
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Figure 3. Example output spectra (histograms) for monoener- 
getic injection at Pqc/Z = 10 18 eV, and fossil jet parameters as 
specified and following particle trajectories as they undergo heli- 
cal motion along field lines. Solid curve shows results for case of 
neglecting gyromotion and simply following the gyrocentre mo- 
tion, and dashed curve shows effect of neglecting time dependence 
of field (valid provided t ma x <t. tdec)- 



3.3 Analytic approach 

No approximations were made in obtaining the energy spec- 
trum, shown as the histogram in Fig. 3, by following particle 
trajectories in the electromagnetic field of a decaying RFP. 
To obtain an approximate analytic result, we can assume 
that a particle follows a helical path with its guiding centre 
on a single field line from one end of the jet to the other. 
Since the induced electric field is proportional to the mag- 
netic field according to Eq. 15, positively charged particles 
will gain energy for pitch angles less than 90° and lose energy 
if their pitch angles are greater than 90° . 

The energy gain on traversing the jet length L will ac- 
tually depend on the pitch angle tp of the helical magnetic 
field line acting as the guiding centre. So, positive particles 
injected into the RFP with r < r cr it will gain energy while 
moving in the positive z direction, and those injected with 
r C rit < r < R will gain energy while moving in the negative 
z direction. The increase in momentum of ultra-relativistic 
particles of charge Ze is 

ZeEL 

Pgain = (20) 



where 



E = 



c costp 



= B { 



tdccCt 



Ji(ar)cj) + Jo(ar)z 

tdccCt 



cosi/> = 
Hence, 

Pgain C 



where 



-^gain 



B 



ZeL B 2 _ ZeB L J (ar) 2 + Ji(ar) 



tdccCt B z ctdccCt Jo(ar) 
o Mar) 2 + Ji(ar) 2 
gain Mar) 



(21) 
(22) 

(23) 
(24) 
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Bo 



10 /ttG I \ Mpc 




Figure 4. (a) Momentum gain of particles injected at one end 
of the RFP of length L and exiting at the other. Solid curve for 
positive particles traveling in positive z direction, dotted for pos- 
itive particles traveling in negative z direction, (b) The function 



R \ I tdec \ 

100 kpc J \GyrJ 



(25) 



and this is plotted in Fig. 4(a). Note that as r — ► r cr it, 

Pgain ^ OO. 

We can work out the energy spectrum as follows, 

(26) 



dN 



dN 



dE s 



dp g 



dN 
dr 



d Pg 



dr 



where dN/dr is the distribution in radius of the injection 
points. For injection at one end of the RFP we would have 
uniform injection over the disk of radius R, giving 

dN 2r „ 

— = -j-r for < r < R. 

dr R 2 

Differentiating Eq. 23 gives 



(27) 



dr 



ZeB L 

tdec 



Ji(ar)Jo(ar) 2 — Ji){ar)J\{ar)M ar ) 



J (ar) 2 



(28) 



which is plotted in Fig. 4(b). 

Thus, from Eqn. 26 we have dN /dp^ iu as a function of 
the parameter r, and from Eqn. 24 we have p ga in as a func- 
tion of the parameter r, and so we can plot dN/dp ga _i n vs. 
Pgain, and this is shown in in Fig. 5. In Fig. 6 we compare the 
analytical spectrum with that obtained by following particle 
orbits. This analytic slope, which asymptotically is E~ 2 , is 
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Figure 5. Spectrum of accelerated particles. (Curves have same 
meaning as in Fig. 4) 

consistent with that of Fig. 5. The agreement reflects a geo- 
metric property of the acceleration mechanism, as particles 
near the reconnection region experience little axial magnetic 
field. By our simplifying assumption that the inductive elec- 
tric field is proportional to the equilibrium magnetic field, 
we find that such particles near the reconnection site then 
are accelerated around the axis by the purely azimuthal elec- 
tric field. This generally keeps them in the structure longer, 
as they do not move radially very much, so they leave only 
when they have drifted down the axis to the end. Those 
which enter the fossil near the reconnection zone will be 
preferentially accelerated, becoming UHECRs. This gives a 
clean mapping of geometry into the energy spectrum. Such 
spectra yield 10 20 eV protons and iron nuclei 26 times more 
energetic. Repeated accelerations in random fossils along a 
particle route could maintain this quality if fossils are com- 
mon enough, and have sufficient decaying energy, in our re- 
gion of the universe. 

Finally, Heinz & Sunyaev (2002) have shown that par- 
ticles should be accelerated at the reverse shock of a micro- 
quasar jet colliding with the interstellar medium, and that 
this may give a contribution to the galactic cosmic rays up to 
about ~10 GeV. Here, we consider the possibility of particle 
acceleration by electric fields in fossil micro-quasars induced 
by the decay of a remnant magnetic bubble after the jets 
have switched off, as in the case of fossil AGN. This mecha- 
nism might work for and pulsar wind nebulae as well as fossil 
micro-quasars, and be responsible for emission in unidenti- 
fied TeV sources, as well as contributing to galactic CR. In 
either case the expected minimum energy would be 




4 DISCUSSION 

Active galaxies have long been considered as a source of 
the UHE CR, for example Rachen & Biermann (1993) con- 
sidered acceleration at shocks associated with hot spots in 




-2-10 1 2 3 
log[(E/Z) / EeV] 



Figure 6. Example output spectra (histograms) for monoener- 
getic injection at poc/Z = 10 18 eV, and fossil jet parameters as 
specified and following particle trajectories as they undergo heli- 
cal motion along field lines. Solid curve shows analytic result. 



lobes of Fanaroff-Riley II (FR II) radio galaxies. There is 
no simple connection between the maximum particle energy 
attainable in an AGN jet, in the hot spot, and in the subse- 
quent bubble. In the case of a jet, either the magnetic field 
through proton synchrotron radiation or pion photoproduc- 
tion interactions will stop the acceleration process (Miicke & 
Protheroe (2001)), whereas in the case of hot spots in lobes 
of FR lis the dimensions of the shock, its speed and magnetic 
field configuration would determine the maximum energy. 
Ultimately, pion photopruduction on the CMBR cannot be 
avoided, and for shock acceleration with plausible parame- 
ters it is difficult to exceed ~10 21 eV (Protheroe (2004)). In 
the case FR II jets (or quasar jets) pion photoproduction on 
the photon field of a luminous accretion disk will cut off the 
spectrum of accelerated protons well below this. However, 
in the case of FR I jets (or the jets of BL Lac objects), hav- 
ing under-luminous disks or advection dominated accretion 
flows, the cut-off will be due to synchrotron losses or pion 
photoproduction on the jet synchrotron photons. The prob- 
lem of escape of the highest energy protons trapped by a 
jet's magnetic field can be solved by the same photoproduc- 
tion interactions, which convert protons to neutrons which 
escape before decaying, and if travelling in the jet direction 
the resulting cosmic ray protons have higher enrgies from 
Doppler boosting (see Protheroe et al. (2003) for the case of 
M87). 

Once accelerated to the ultra-high energy region, pro- 
tons can more readily escape their birth fossil. They will 
move through the tens of Mpc where the GZK photo-pion 
scattering from the background microwave radiation begins 
to take its toll. As well, the UHECRs will scatter from other 
fossils. This can make them appear isotropic if their direc- 
tions are scattered enough, and they can be further acceler- 
ated by these other fossils as well. 

Our model is geometrically simple, to capture the essen- 
tial geometry of reverse field pinches, and the acceleration 
physics that affords to particles, depending on where they 
are in the structure. We have neglected important effects: 

(i) Scattering - this in principle will rearrange particle 
orbits, moving some out of the high-acceleration zone by 
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scattering their pitch angles, but also scattering others into 
the preferential region. In the highly ordered field of AGN 
fossils, this may be less important ant for energization to 
UHE of a pre-existing population of high energy cosmic rays 
accelerated, e.g. in the host galaxy of the dead AGN. Fu- 
ture work should nevertheless include investigation of scat- 
tering of lower energy particles by the magnetic turbulence 
expected in such sources, generated by the inductive effects 
of their slow collapse and by the ongoing reconnection events 
near the axial field reversal sites. 

(ii) Energy losses - here we have neglected energy losses 
by, e.g. pion photoproduction and adiabatic expansion. For 
the range of model parameters used in the present work, 
these occur more slowly than acceleration, but should be 
included in future work. 

(iii) Geometry. Assuming a cylindrical reversed field 
pinch omits the many possible curvatures a more spherical 
RFP i.e., a spheromak would entail. External pressure will 
not be isotropic, generally. But if a fossil becomes detached 
from its parent galaxy presumably by reconnection of fields 
at the feet of the former jet structure it can rise, buoyant, 
in any cluster gas. The shape will resemble a fat torus, with 
the currents forming loops. We simplified this, feeling that 
torodial curvature would complicate our calculations with- 
out adding new physics. This can be relaxed in future, using 
more complicated equilibria. 

(iv) Initial energy. To isolate the highest energy particle 
issues, we began with a population of constant energy par- 
ticles, e.g. 1 EeV, but provided the injection energies are 
much less than the minimum energy gained, by particles be- 
ing accelerated along the fossil jet, axis the final spectral 
shape is insensitive to the initial energy. Generalizing to an 
initial energy spectrum will be done in future. This said, we 
have produced spectra that yield very high energy particles. 
10 20 eV protons and iron nuclei 26 times more energetic 
still. If 1% of the universe is taken up with fossils (Kron- 
berg (1994)), of size < 1 Mpc, then there should be many 
fossils within 100 Mpc of us. A fresh structure is under con- 
struction at M87, 16 Mpc away. This leads us to predict 
that while Auger should see a smooth UHECR background 
above the GZK cutoff, ~6 x 10 19 cV. Au ger may see clumps 
of UHECR in the sky, which point back to nothing specially 
luminous. These will be fossils we didn't know about-plus 
some we did. The UHECR will then become a diagnostic of 
regions storing vast magnetic energies, yet perceptible only 
through equally energetic particles, bringing word of them 
across great distances. 

The total energy stored in the magnetic field of an RFP 



U B = 1.3 x 10 59 ( i 



10 jiiG 



R 



( * 

^ Mpc J \ 100 kpc 



(erg). (30) 



These magnetic structures must have survived since the 
quasar was last active, implying a decay time scale idee of 
the order of Gyr, so 



Ub = 

Cdec 



U B = (4.0 x 10 42 erg s" 1 ) 



Bn 



10 HG 



L 



R 



1 Mpc I \ 100 kpc 



^dec \ 

GyTj 



(32) 



assuming that the reconnection is continuous and uniform 
along the cylindrical reconnection region, such that the 
whole magnetic structure slowly decays exponentially on 
time scale idee while maintaining the stable configuration 
given by Eq. 3. 

O'Neill et al. (2006) have examined the interactions of 
AGN jets with magnetized cluster media using MHD jet 
simulations in a cluster-like atmosphere extending over dis- 
tances ~200 kpc. They find that steady jets asymptotically 
deposit ~50% of their power into ambient thermal energy, 
~25% into heating of the ICM, and the remaining energy to 
be mostly stored in the jet backflow cocoon, and that mag- 
netic energy introduced by the jet is generally amplified, 
especially by sheared flows; see also Ryu et al. (1998). 

For an extragalactic source distribution locally produc- 
ing an E~ 2 spectrum of protons, Lipari (2005) estimates the 
local power requirement to be ~ 10 50 erg Mpc -3 y _1 . De- 
caying magnetic fields with local filling factor tjb lose energy 
at a rate 



ub 



10 53 



'rjB 



Bp \ ( tdec \ 

10 »G) \Gyr) 



erg Mpc 3 y 1 



Magnetic fields from quasars can fill up to 5-20% of the 
intergalactic medium (Furlanetto & Loeb (2001)) - proba- 
bly higher locally since our Galaxy is in a "Wall". Indeed, 
Copal-Krishna & Wiita (2001) estimate the fractional rel- 
evant volume that radio lobes born during the quasar era 
cumulatively cover is ~0.5. Hence our crude energetics ar- 
guments show fossil AGN structure decay could well be re- 
sponsible for the observed UHE CR. 

The spectrum of accelerated particles will cut off at 
some maximum momentum determined by either the finite 
thickness of the reconnection zone (recall that in the an- 
alytic approximation as r -» r cr it, Pgain — » oo), or by the 
gyroradius increasing so that it is no longer much less than 
the radius of the fossil. From Fig. 1, we see that for r < R 
the magnetic field is in the range 0.4£?o < B < Bo. Hence, 
the condition tl <C R implies 



< {10 21 Z eV 



R 



Bo 



10 juG I \ 100 kpc 



(33) 



(31) 



The UHECR spectrum shows a steepening just below 
10 18 eV (the "dip") followed by a flattening at around 10 19 
eV (the "ankle") before the possible observed GZK cut-off 
just above 10 20 eV. Berezinsky et al. (2006) favour a scenario 
in which extragalactic cosmic rays are injected with an E~ 2 ' 7 
spectrum above above ~2 x 10 ls eV, and show that the ob- 
served spectal feature referred to as the "dip" could then 
be due to to Bethe-Heitler pair production on the CMBR. 
While their model fits are impressive, we believe that mod- 
els in which the galactic cosmic rays extend up to a few 
10 18 eV and have extragalacic cosmic rays injected with an 
E~ spectrum, or similar, e.g. following diffusive shock ac- 
celeration at Fanaroff-Riley Class II radio galaxy hot-spots 
(Rachen & Biermann 1993), are by no means ruled out. In- 
deed, Hillas (2006) has shown that by adding galactic and 
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extragalactic spectra with various power-law source spectra 
and compositions several possibilities, including source spec- 
tra with an E~ 2 ' 2 spectrum, are consistent with the data. 
While a single ideal cylindrical RFP with particle injection 
over its ends yields an E~ 2 spectrum, the spectrum due to 
several with different maximum energies would inevitably 
be steeper. Very recent results from the Pierre Auger Ob- 
servatory (Yamamoto et al. 2007) are not yet conclusive, 
but appear to be consistent with source spectral index in 
the range 2.2-2.5. 

The spectrum of UHECR observed at Earth would have 
contributions from nearby fossil jets at different distances, 
with different powers and each having different dimensions 
and magnetic fields, and hence a range of E^ in and E™^*. 
For an individual fossil, the cut-off is expected to be rigidity 
dependent, implying the observed composition would change 
from light to heavy close to the cut-off if one or two nearby 
AGN fossils dominate. However, if distant sources domi- 
nate nuclei will be photo-disintegrated by interactions with 
CMBR photons, and in this case the composition would re- 
main light to the highest energies if distant sources or fossils 
dominated. Otherwise the composition could be mixed near 
the observed cut off. 

Nearby potential sources of UHE CR would include the 
bubbles around Cen A, the nearest powerful radio galaxy, 
and around M87. Radio observations at 20 cm of Cen A by 
Morganti et al. (1999) indicate the presence of large features, 
one of which has been interpreted by Saxton et al. (2001) 
as a buoyant bubble of plasma deposited by an intermit- 
tently active jet. Even closer, in M87 there are two bubbles 
of synchrotron emission which appear to be inflated by out- 
flows from the central engine visible in the low frequency 
radio data of Owen et al. (2000), with an amazingly sharp 
edge indicating the possible presence of a shock. Given the 
very high mass of the central black hole in M87, its jet, al- 
though fairly powerful, is weak compared to its Eddington 
limit, making M87's bubbles a prime cadidate as a UHE CR 
accelerator. 

As we expect most of the fossil jets to be below the 
sensitivity of current radio telescopes, it is impossible at 
the present time to make firm predictions for the expected 
UHECR intensity at Earth. A low synchrotron surface 
brightness, only observable if at all at low frequency, does 
not necessarily mean that the magnetic field is low. Rather, 
we would argue that it is the cooling of the relativistic elec- 
trons by synchrotron emission, since the last AGN outburst 
which filled the bubble, that has reduced the synchrotron 
cmissivity. However, such fossils may well observable when 
the SKA (http://www.skatelescope.org/) is commissioned, 
hopefully within ~10 years. Nevertheless, we have demon- 
strated that it is possible for this process to accelerate 
protons to UHE, and nuclei to a Z times higher energy, 
and shown that the power requirements may reasonably be 
achieved given plausible volume filling factors. 

Very recent anisotropy data from the Pierre Auger Ob- 
servatory (Armengaud et al. (2007) ,Harari et al. (2007)) 
show no significant large-scale anisotropy or correlation of 
cosmic ray arrival directions with BL Lac objects. Perhaps 
this is not surprising because propagation of cosmic rays 
through intergalactic magnetic fields will tend to isotropize 
their arrival directions, except perhaps for UHE CR pro- 
tons at the highest energies from very nearby extragalactic 



sources. With longer exposure times Auger may see such 
anisotropics, in which case the most probable source candi- 
dates would be bubbles around Cen A and M87. Neurino ob- 
servatories may, in principle, see UHE neutrinos originating 
from directions close to UHE cosmic ray acceleration sites, 
through pion photoproduction interactions of the UHE CR 
on the CMBR during propagation from their sources. These 
"GZK neutrinos" would appear as haloes around the true 
source direction. Of course, if such haloes appear where there 
are no obvious source candidates, they could well be from 
fossil jets with very low surface radio brightness. 

Future work will include exploring observational conse- 
quences of this model at radio, X-ray and gamma-ray en- 
ergies. Will electron acceleration (in opposite direction to 
protons) be killed by streaming instability, or by other pro- 
cesses? If so, this will heat the plasma. Could this be seen 
in X-rays? If not, radio synchrotron and gamma-ray sig- 
nals could be generated. Could these be detectable by ex- 
isting radio telescopes or by the SKA in coming decades, or 
by GLAST (http://glast.stanford.edu/)? How much of the 
energy of a fossil AGN goes into cosmic rays? Finally, we 
mention the possibility that this mechanism operating for 
microquasars might be responsible for acceleration of some 
of the Galactic population of cosmic rays. 

In conclusion, remnants of jets and their surrounding 
cocoons may still be present around or close to galaxies 
which contain AGN which are now no longer active. These 
fossil jets are colossal MHD structures and may have total 
energies ~ 10 60 erg. After a jet switches off, the fossil jet 
will organize itself into a stable configuration, such as the 
reverse-field pinch in which the longitudinal magnetic field 
changes direction at a critical distance from the axis, lead- 
ing to re-connection there, and to slow decay of the large- 
scale RFP field. We have shown that this field decay in- 
duces large-scale electric fields which accelerate cosmic rays 
an E~ 2 power-law up to ultra-high energies. Energetics ar- 
guments show that this provides a plausible mechanism for 
the origin of the UHE CR. 
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